Color  xc   yc  (LatLon LCarr)  Area Length Branches (  x0/y0     x1/y1  )  Tilt Location
      1 1204 1712 (N40E06 L= 21)   505     87        1 (1166/1685-1235/1740)  45.0 E08-E03 N37-N43
      2 1296 1633 (N37W02 L= 30)   311     58        0 (1272/1605-1319/1660) -45.0 W00-W04 N34-N38
      3  425 1551 (N12E53 L=335)  1172    240        1 ( 400/1461- 445/1645) -63.4 E58-E48 N06-N18
      4  959 1423 (N14E13 L= 15)   511    156        0 ( 936/1336- 992/1493) -52.1 E15-E08 N09-N18
      5 1736 1365 (N32W40 L= 68)   365    133        0 (1708/1329-1774/1396)  49.4 W37-W44 N28-N35
      6  308 1263 (S07E52 L=335)   294     62        0 ( 287/1233- 332/1298)  63.4 E53-E51 S08-S04
      7 1752 1223 (N23W43 L= 70)   368     92        0 (1729/1172-1775/1270) -45.0 W40-W45 N20-N25
      8  523 1102 (S13E32 L=356)   784    203        1 ( 460/1008- 583/1195)  69.0 E34-E29 S19-S06
      9  142 1114 (S17E66 L=321)   674    151        2 ( 126/1034- 157/1168) -48.8 E70-E63 S21-S13
     10 1134  972 (S08W07 L= 35)   931    269        2 (1038/ 843-1246/1104)  74.5 W04-W10 S16-N02
     11 1677  737 (S07W48 L= 76)  1221    268        1 (1573/ 659-1761/ 860)  61.9 W42-W51 S13-N02
     12  845  806 (S23E06 L= 22)   642    142        1 ( 774/ 784- 921/ 832)  24.0 E10-E01 S25-S20
  

    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