Color  xc   yc  (LatLon LCarr)  Area Length Branches (  x0/y0     x1/y1  )  Tilt Location
      1  659 1565 (N21E38 L=181)   857    274        0 ( 616/1486- 701/1665) -32.7 E45-E31 N16-N26
      2  513 1595 (N20E50 L=168)   444    109        1 ( 464/1575- 549/1624)  38.7 E53-E47 N18-N22
      3 1600 1259 (N25W32 L=250)   464     92        2 (1576/1212-1626/1300) -68.2 W29-W34 N22-N27
      4 1231 1140 (N09W09 L=227)   318     49        1 (1201/1114-1254/1163)   0.0 W07-W10 N07-N10
      5  158 1013 (S21E63 L=155)   778    177        2 ( 128/ 973- 195/1058)  45.0 E66-E59 S24-S18
      6  347  987 (S19E44 L=174)   351    139        0 ( 331/ 926- 379/1040)  90.0 E45-E41 S22-S15
      7 1824  990 (N15W56 L=274)   792    200        0 (1771/ 957-1872/1016)   9.5 W50-W62 N11-N16
      8 1736  958 (N10W48 L=266)   431    111        0 (1707/ 924-1767/1006) -33.7 W44-W51 N08-N12
      9  620  852 (S22E21 L=197)   667    302        2 ( 552/ 746- 696/ 978)  70.6 E24-E18 S29-S12
     10 1177  819 (S11W14 L=232)   417    157        2 (1125/ 771-1234/ 876)  60.3 W10-W16 S14-S07
     11 1665  545 (S15W57 L=276)   991    117        2 (1625/ 511-1710/ 578)  39.8 W52-W60 S16-S11
     12 1046  501 (S33W16 L=234)   429     73        3 (1014/ 476-1085/ 528)  36.9 W13-W18 S35-S31
     13 1512  365 (S29W55 L=273)   685     86        1 (1485/ 339-1541/ 405) -59.0 W51-W57 S31-S26
  

    OPENR, u, binfile, /COMP, /GET_LUN
    head = FLTARR(6) ;[xsize,ysize,radius,hh,mm,ss]
    READU, u, head
    data = INTARR(FIX(head[0]),FIX(head[1]))
    READU, u, data
    FREE_LUN, u
    

• each filament has the value of Fil_ID
• all values are in pixel

  c = WHERE(data EQ Fil_ID, area)
  f = '(2(A,I6))'
  PRINT, FORM=f, "Area of Fil_ID", Fil_ID," =", area 

  data = ROTATE(data,7)
  Solar_y = 2048-Solar_y