Color  xc   yc  (LatLon LCarr)  Area Length Branches (  x0/y0     x1/y1  )  Tilt Location
      1  655 1605 (N20E41 L=180)   493    107        2 ( 611/1561- 687/1632) -12.5 E45-E36 N17-N21
      2 1330 1603 (N40W03 L=223)   450     89        1 (1282/1588-1381/1619)  18.4 W00-W06 N37-N41
      3  914 1556 (N23E22 L=198)   370    200        0 ( 893/1483- 972/1611) -26.6 E25-E15 N20-N27
      4 1253 1428 (N25W03 L=223)   413    118        0 (1221/1373-1291/1471)  82.9 W01-W04 N20-N28
      5  673 1129 (S07E23 L=198)   201     73        1 ( 640/1098- 703/1156) -11.3 E25-E20 S08-S06
      6  213 1058 (S22E59 L=161)  1049    136        3 ( 190/ 998- 229/1128) -76.0 E62-E56 S25-S17
      7  148  955 (S29E69 L=152)  1605    249        3 ( 115/ 828- 179/1035)  39.5 E83-E63 S37-S23
      8  443  900 (S26E34 L=186)   467    115        0 ( 395/ 869- 497/ 938)  -6.3 E38-E29 S27-S24
      9  725  778 (S26E11 L=209)   543    209        0 ( 683/ 722- 787/ 866) -36.9 E15-E06 S28-S21
     10 1265  717 (S15W22 L=243)   279     54        2 (1243/ 687-1288/ 744) -56.3 W20-W23 S16-S13
     11  311  575 (S52E52 L=169)   513     68        0 ( 295/ 546- 329/ 608) -76.0 E55-E49 S53-S49
     12 1404  454 (S26W43 L=264)   344     77        0 (1376/ 422-1437/ 492)  90.0 W42-W44 S27-S22
  

    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