Color  xc   yc  (LatLon LCarr)  Area Length Branches (  x0/y0     x1/y1  )  Tilt Location
      1  301 1455 (N05E60 L=132)  2601    502       10 ( 207/1376- 384/1571)  59.7 E71-E53 N00-N13
      2  865 1377 (N10E17 L=175)   426    140        1 ( 835/1303- 890/1432) -63.4 E19-E13 N05-N13
      3  116 1157 (S14E72 L=120)  2835    291        5 (  88/1088- 147/1233) -58.0 E77-E67 S17-S09
      4  227 1178 (S12E58 L=134)   356     82        1 ( 204/1142- 240/1205) -76.0 E60-E56 S14-S10
      5  321  932 (S26E45 L=147)  2241    507        6 ( 218/ 774- 440/1104)  55.2 E54-E37 S36-S13
      6 1267  897 (S09W17 L=209)   409     96        1 (1213/ 874-1309/ 923)   9.5 W13-W19 S10-S07
      7 1276  754 (S16W21 L=213)  1170    337        2 (1132/ 656-1389/ 874)  63.4 W14-W25 S25-S07
      8  917  752 (S25E00 L=192)   622    120        1 ( 847/ 704- 973/ 805)  63.4 E03-W02 S29-S21
      9  326  618 (S46E46 L=146)   372     79        1 ( 298/ 599- 359/ 648)  -9.5 E48-E42 S47-S43
     10 1567  476 (S24W52 L=244)  3155    742        2 (1275/ 366-1845/ 580)  35.9 W36-W76 S38-S09
     11  550  545 (S47E24 L=168)   352     45        2 ( 527/ 522- 577/ 569) -21.8 E26-E21 S48-S44
     12 1321  512 (S29W33 L=225)   332     64        1 (1283/ 497-1361/ 526)  36.9 W30-W35 S30-S27
  

    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