Color  xc   yc  (LatLon LCarr)  Area Length Branches (  x0/y0     x1/y1  )  Tilt Location
      1  976 1743 (N35E27 L=299)   342     50        0 ( 961/1708- 993/1775) -59.0 E29-E24 N32-N37
      2  859 1697 (N28E31 L=295)   461    164        1 ( 785/1649- 914/1723)  -4.8 E37-E25 N26-N29
      3 1247 1439 (N22W02 L=328)   279    115        0 (1194/1407-1290/1460)  -9.5 E02-W04 N21-N23
      4  851 1236 (N00E15 L=311)   388    140        1 ( 810/1184- 897/1291) -23.2 E18-E11 S01-N02
      5 1634 1076 (N13W35 L=  1)   559    109        1 (1615/1024-1652/1131) -81.9 W33-W36 N09-N16
      6  661 1096 (S12E23 L=303)   254    124        0 ( 607/1075- 717/1111)  14.0 E26-E18 S13-S11
      7 1784  977 (N13W50 L= 16)   816    220        1 (1748/ 881-1825/1078) -36.9 W44-W56 N08-N17
      8  847  868 (S20E06 L=320)   283    104        0 ( 795/ 855- 899/ 881)  26.6 E09-E03 S21-S18
      9  601  835 (S28E21 L=305)   346    120        0 ( 540/ 824- 663/ 848)  18.4 E25-E16 S30-S26
     10  787  648 (S35E04 L=322)   867    205        5 ( 687/ 614- 881/ 689)  36.0 E10-W01 S38-S30
     11 1481  397 (S27W54 L= 20)  1281    172        2 (1396/ 363-1555/ 439)  41.6 W48-W58 S31-S23
     12 1249  254 (S44W51 L= 17)   524     93        0 (1223/ 234-1290/ 281)  76.0 W46-W54 S45-S41
  

    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