Color  xc   yc  (LatLon LCarr)  Area Length Branches (  x0/y0     x1/y1  )  Tilt Location
      1 1392 1796 (N39W59 L=287)  2250    243        4 (1284/1759-1516/1820) -12.5 W46-W73 N34-N42
      2 1655 1610 (N20W64 L=292)  2590    302        3 (1557/1521-1779/1719) -70.6 W60-W72 N11-N28
      3 1278 1613 (N31W34 L=262)   605    180        0 (1194/1593-1364/1633) -28.6 W28-W39 N27-N33
      4  984 1445 (N29W08 L=236)   204     75        0 ( 948/1423-1020/1464) -45.0 W06-W10 N26-N30
      5  598 1413 (N37E19 L=209)   889    346        0 ( 464/1367- 780/1457)  -6.6 E30-E04 N34-N38
      6 1258 1281 (N12W20 L=247)   496    163        0 (1218/1208-1302/1349) -90.0 W18-W21 N07-N17
      7 1766 1104 (S10W50 L=277)   333     73        0 (1739/1071-1786/1143) -78.7 W48-W51 S12-S07
      8 1836 1089 (S13W57 L=284)   375     83        0 (1801/1062-1870/1118)   6.3 W52-W61 S14-S12
      9  200  893 (N14E61 L=166)   858     90        0 ( 155/ 861- 245/ 925) -21.8 E67-E56 N11-N16
     10  869  812 (S05E14 L=214)   798    329        1 ( 764/ 692- 929/ 867) -71.6 E19-E10 S12-N00
     11 1828  765 (S33W53 L=281)   790    124        1 (1807/ 722-1864/ 827) -42.0 W50-W60 S37-S28
     12 1680  702 (S32W37 L=264)   647    124        0 (1633/ 666-1721/ 744)  14.0 W32-W40 S34-S29
     13 1166  627 (S22E01 L=226)   952    240        2 (1094/ 533-1205/ 723)  39.3 E08-W03 S26-S17
     14 1573  515 (S42W26 L=254)  1338    185        2 (1495/ 459-1635/ 591)  19.7 W18-W32 S45-S38
     15  531  375 (S23E54 L=174)   477     90        0 ( 517/ 330- 549/ 424)  36.9 E57-E48 S26-S20
     16  649  305 (S31E50 L=178)   598    133        0 ( 605/ 260- 677/ 362)  19.7 E56-E42 S33-S27
     17 1140  246 (S50E18 L=209)  1710    281        3 (1007/ 194-1266/ 299)  -2.1 E32-E05 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