Color  xc   yc  (LatLon LCarr)  Area Length Branches (  x0/y0     x1/y1  )  Tilt Location
      1  764 1706 (N30E34 L= 66)  1135    239        2 ( 663/1640- 848/1779) -15.9 E46-E25 N26-N34
      2  437 1681 (N24E61 L= 39)  1584    154        3 ( 399/1607- 494/1723)  90.0 E68-E57 N18-N27
      3 1137 1649 (N34E07 L= 93)   691    107        1 (1103/1603-1171/1687) -50.2 E10-E03 N30-N37
      4 1120 1267 (N09W01 L=101)  1197    529        1 (1076/1066-1206/1485) -54.5 E06-W09 S01-N21
      5 1321 1399 (N21W10 L=111)   245    100        0 (1281/1377-1373/1423)   0.0 W07-W13 N20-N22
      6 1784 1287 (N27W48 L=148)   854    118        0 (1757/1234-1809/1348) -66.8 W45-W50 N23-N30
      7 1335  813 (S12W23 L=123)   371     83        0 (1293/ 776-1365/ 852) -18.4 W19-W25 S14-S10
      8 1457  781 (S11W31 L=132)   548    184        0 (1393/ 746-1545/ 816)  33.7 W27-W36 S14-S07
      9  974  722 (S25W03 L=104)   375    117        2 ( 929/ 680-1016/ 758) -18.4 W00-W06 S27-S23
     10 1716  647 (S12W56 L=156)   345     57        1 (1697/ 626-1737/ 670) -71.6 W53-W58 S13-S10
     11  899  598 (S35W01 L=102)   214     61        0 ( 867/ 583- 931/ 619) -14.0 E01-W03 S35-S33
     12 1492  414 (S31W51 L=152)   421     92        0 (1448/ 397-1538/ 433)   5.7 W46-W56 S31-S29
     13  846  196 (S68W16 L=117)   861    122        2 ( 788/ 162- 905/ 231)  -2.5 W06-W29 S70-S65
  

    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