Color  xc   yc  (LatLon LCarr)  Area Length Branches (  x0/y0     x1/y1  )  Tilt Location
      1 1502 1719 (N31W59 L=232)  1145    198        2 (1397/1694-1579/1742) -19.4 W49-W67 N27-N34
      2 1130 1665 (N39W27 L=199)   986    243        0 (1016/1618-1230/1702)  -6.0 W16-W35 N37-N39
      3 1336 1446 (N20W30 L=202)   270     85        1 (1297/1423-1370/1470) -68.2 W28-W31 N17-N22
      4  365 1287 (N33E40 L=132)   922    171        2 ( 278/1265- 446/1305) -15.5 E50-E32 N30-N36
      5  907 1184 (N15E03 L=169)   367    121        1 ( 847/1160- 963/1203) -39.8 E06-W00 N12-N17
      6  527 1161 (N22E28 L=145)   444    119        1 ( 483/1118- 584/1181)   7.1 E31-E23 N20-N22
      7 1776 1032 (S14W49 L=222)  1257    211        5 (1713/ 955-1847/1125)  23.6 W42-W58 S18-S10
      8 1664 1016 (S12W39 L=212)   562    236        1 (1593/ 945-1708/1117)  30.3 W32-W44 S15-S07
      9  181  972 (N17E62 L=111)  2010    249        3 (  87/ 918- 255/1045)   5.9 E81-E52 N14-N20
     10 1111  950 (S03W03 L=176)   242     98        0 (1060/ 939-1159/ 963) -31.0 W00-W06 S03-N00
     11 1841  893 (S25W55 L=228)  1081     96        2 (1804/ 856-1878/ 940)  11.3 W50-W60 S26-S21
     12 1407  736 (S22W16 L=189)   275     61        1 (1393/ 696-1423/ 771)  63.4 W14-W17 S24-S20
     13 1210  627 (S24W02 L=174)   790    250        1 (1154/ 523-1270/ 742)  42.3 E04-W07 S29-S18
     14 1015  609 (S20E10 L=162)   869    258        4 ( 976/ 488-1064/ 722)  47.3 E16-E04 S27-S14
     15  697  479 (S21E35 L=138)  1461    395        2 ( 646/ 307- 759/ 647)  38.4 E48-E24 S31-S11
     16 1271  540 (S31W04 L=176)   384    111        1 (1217/ 509-1318/ 577)   8.1 W00-W07 S32-S29
     17 1520  404 (S50W22 L=195)  1083    178        0 (1432/ 373-1592/ 442)  -3.2 W12-W30 S52-S47
  

    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