Color  xc   yc  (LatLon LCarr)  Area Length Branches (  x0/y0     x1/y1  )  Tilt Location
      1  530 1677 (N58E27 L=259)  2067    337        4 ( 374/1591- 650/1746)   7.9 E49-E11 N54-N61
      2 1306 1554 (N27W32 L=318)   291     40        0 (1282/1537-1334/1575) -56.3 W30-W34 N25-N28
      3 1385 1475 (N20W34 L=320)   445     95        2 (1337/1444-1418/1506) -18.4 W30-W36 N18-N22
      4  802 1428 (N32E04 L=282)   251     54        0 ( 785/1400- 825/1458) -76.0 E05-E03 N29-N33
      5  242 1348 (N39E57 L=229)   428     89        1 ( 208/1329- 281/1366)  -5.2 E62-E51 N38-N41
      6 1157 1213 (N11W12 L=298)   229     62        0 (1128/1194-1185/1237)  14.0 W09-W13 N09-N11
      7 1663 1016 (S12W39 L=325)   489     78        0 (1649/ 959-1682/1061)  71.6 W37-W41 S15-S09
      8  229  759 (N04E62 L=224)   642     85        0 ( 211/ 715- 250/ 803)  59.0 E63-E58 N01-N06
      9 1616  681 (S31W32 L=318)   384    103        0 (1562/ 664-1662/ 695) -20.6 W27-W35 S33-S30
     10  499  610 (S10E43 L=243)   459    107        2 ( 472/ 561- 525/ 655)  35.5 E46-E39 S12-S07
     11 1007  404 (S34E18 L=268)  1293    386        0 ( 860/ 322-1133/ 534)  13.5 E30-E05 S37-S28
     12 1207  312 (S47E08 L=278)   296     45        0 (1179/ 298-1236/ 330)   9.5 E10-E04 S47-S45
  

    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