Color  xc   yc  (LatLon LCarr)  Area Length Branches (  x0/y0     x1/y1  )  Tilt Location
      1  591 1574 (N24E43 L=168)  2931    515        5 ( 446/1462- 801/1730) -12.9 E72-E23 N19-N31
      2  958 1439 (N22E13 L=198)   244     83        0 ( 925/1421-1005/1458)   0.0 E15-E09 N21-N23
      3 1466 1358 (N28W23 L=233)   289    121        0 (1429/1312-1502/1397) -33.7 W19-W25 N25-N29
      4  760 1310 (N11E22 L=189)   342    122        0 ( 709/1270- 807/1350) -15.9 E25-E18 N08-N12
      5  727 1215 (N04E21 L=189)   866    357        4 ( 601/1171- 908/1265)  27.9 E28-E11 N00-N09
      6 1102 1120 (N06W03 L=213)   352    196        1 (1076/1036-1128/1213)  90.0 W01-W03 N00-N11
      7  504  893 (S19E31 L=180)   548    160        1 ( 480/ 818- 534/ 960)  90.0 E32-E28 S23-S14
      8  601  757 (S25E22 L=188)   386    113        2 ( 549/ 736- 648/ 778)   7.1 E26-E18 S26-S23
      9 1820  674 (S06W69 L=279)   492     35        1 (1799/ 658-1840/ 695)  -8.1 W65-W72 S07-S04
     10 1656  473 (S21W62 L=273)   472     66        0 (1627/ 450-1691/ 498)  45.0 W59-W64 S23-S19
     11 1278  371 (S37W36 L=247)   328     57        1 (1255/ 344-1301/ 399) -26.6 W33-W39 S38-S34
     12 1418  376 (S33W47 L=258)   495     82        2 (1377/ 360-1456/ 398)  15.9 W43-W50 S34-S31
  

    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