Color  xc   yc  (LatLon LCarr)  Area Length Branches (  x0/y0     x1/y1  )  Tilt Location
      1 1410 1712 (N52W05 L= 79)   539     72        0 (1376/1677-1439/1742)  71.6 W02-W07 N47-N56
      2 1144 1454 (N23E05 L= 69)   797    334        1 (1098/1316-1195/1563) -73.6 E08-E03 N13-N30
      3  766 1507 (N16E29 L= 45)   267     36        1 ( 745/1480- 784/1527) -26.6 E31-E27 N14-N17
      4 1741 1385 (N39W42 L=116)   771     97        1 (1699/1351-1775/1407)   0.0 W37-W46 N37-N41
      5 1302 1293 (N18W09 L= 82)   315    115        1 (1253/1244-1340/1317)  -9.5 W05-W11 N15-N19
      6 1518 1131 (N15W26 L=100)   252     32        0 (1502/1108-1535/1161) -71.6 W24-W27 N13-N16
      7  166 1104 (S21E67 L=  7)   984     98        2 ( 147/1056- 189/1151) -71.6 E69-E63 S23-S17
      8  239 1089 (S20E56 L= 17)   974    136        3 ( 209/1023- 284/1151)  90.0 E59-E51 S24-S16
      9 1629 1056 (N14W36 L=109)   339     68        0 (1609/1028-1645/1084) -26.6 W33-W37 N12-N15
     10  747  641 (S34E06 L= 68)   217     38        0 ( 727/ 626- 777/ 659)  14.0 E07-E03 S35-S33
     11  719  351 (S57W02 L= 76)   307     57        0 ( 693/ 322- 747/ 377) -15.9 E01-W06 S58-S54
     12 1058  222 (S53W41 L=115)  1771    237        3 ( 979/ 168-1152/ 285)  -4.2 W28-W57 S55-S49
  

    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