Color  xc   yc  (LatLon LCarr)  Area Length Branches (  x0/y0     x1/y1  )  Tilt Location
      1  916 1903 (N73W47 L= 76)   915     93        0 ( 894/1867- 941/1944)   8.1 W33-W83 N68-N75
      2 1258 1506 (N30W28 L= 58)   308    115        0 (1211/1480-1315/1521) -38.7 W25-W31 N26-N31
      3 1739 1500 (N15W67 L= 97)   486     68        2 (1703/1482-1772/1514) -11.3 W62-W72 N13-N15
      4 1131 1422 (N28W17 L= 46)   261     69        1 (1106/1386-1150/1461)  78.7 W14-W18 N25-N30
      5 1686 1421 (N13W56 L= 86)   426    110        1 (1648/1370-1710/1461)  81.9 W54-W59 N09-N16
      6 1747 1423 (N11W63 L= 92)   515     52        1 (1722/1399-1773/1450) -76.0 W59-W64 N08-N12
      7 1074 1314 (N23W10 L= 39)   504    190        1 (1018/1231-1111/1384)  33.7 W04-W13 N19-N26
      8  705 1355 (N34E14 L= 16)   237     72        1 ( 672/1324- 733/1378) -53.1 E15-E11 N31-N35
      9 1756 1228 (N00W54 L= 83)   881    128        0 (1715/1169-1793/1285) -74.1 W51-W56 S03-N04
     10  774 1148 (N19E13 L= 17)   237     53        0 ( 750/1127- 804/1173) -45.0 E13-E10 N18-N21
     11  614 1071 (N18E24 L=  5)   250     98        0 ( 577/1033- 651/1095) -51.3 E26-E22 N15-N20
     12 1654  888 (S16W37 L= 67)   567    122        1 (1599/ 859-1702/ 934)   0.0 W32-W41 S17-S14
     13 1535  668 (S27W25 L= 55)   619     95        2 (1511/ 625-1559/ 720)  80.5 W23-W26 S29-S23
     14  645  663 (S07E31 L=359)   234     79        0 ( 628/ 627- 678/ 693)  36.9 E32-E28 S08-S05
     15  598  450 (S19E43 L=347)   413     79        0 ( 558/ 432- 641/ 464)  -8.1 E46-E39 S20-S18
     16  879  421 (S27E24 L=  5)   697    154        1 ( 813/ 389- 953/ 457)   4.8 E30-E18 S28-S25
  

    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