Color  xc   yc  (LatLon LCarr)  Area Length Branches (  x0/y0     x1/y1  )  Tilt Location
      1 1343 1624 (N31W36 L=138)   335     57        1 (1316/1594-1371/1652) -76.0 W34-W37 N29-N33
      2 1628 1536 (N20W54 L=156)   461     28        1 (1606/1512-1651/1561) -56.3 W52-W56 N18-N22
      3  553 1390 (N32E25 L= 77)   363    127        0 ( 501/1338- 604/1422)  26.6 E29-E21 N29-N33
      4 1815 1022 (S13W55 L=157)   319     39        0 (1797/ 993-1831/1044) -56.3 W52-W56 S14-S11
      5 1518  993 (S09W29 L=132)   415     82        1 (1496/ 952-1537/1039)  51.3 W27-W31 S11-S06
      6  297  864 (N04E51 L= 51)  1033    163        0 ( 231/ 828- 348/ 928)  17.1 E58-E45 N02-N07
      7 1064  831 (S10E01 L=101)   799    221        2 (1005/ 728-1122/ 913)  41.6 E06-W03 S15-S06
      8  478  749 (S05E39 L= 63)   499    115        2 ( 449/ 694- 522/ 777)  29.7 E42-E35 S08-S03
      9 1815  699 (S34W57 L=159)   995     85        1 (1787/ 652-1844/ 744) -41.2 W52-W63 S37-S30
     10 1705  664 (S34W44 L=146)   345     82        0 (1677/ 632-1731/ 695)  31.0 W40-W46 S35-S31
     11  494  618 (S13E42 L= 60)   827    109        2 ( 462/ 582- 533/ 670) -78.7 E45-E39 S15-S10
  

    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