Color  xc   yc  (LatLon LCarr)  Area Length Branches (  x0/y0     x1/y1  )  Tilt Location
      1  727 1734 (N32E42 L=351)  1540    301        1 ( 682/1659- 785/1827) -36.3 E50-E34 N28-N39
      2  846 1745 (N36E35 L=359)   399    111        0 ( 818/1710- 880/1780) -26.6 E39-E31 N34-N38
      3  694 1542 (N19E35 L=359)  1806    413        5 ( 575/1387- 749/1625) -31.4 E45-E27 N09-N23
      4 1148 1538 (N30E06 L= 28)   745    196        1 (1058/1464-1233/1605) -15.9 E13-W01 N27-N33
      5  496 1470 (N10E45 L=349)   968    294        2 ( 466/1381- 554/1576) -45.0 E50-E37 N05-N16
      6 1626 1428 (N37W31 L= 65)   613    116        1 (1576/1379-1684/1460) -11.3 W26-W36 N35-N38
      7 1342 1369 (N25W11 L= 44)   527    142        1 (1281/1324-1427/1408)   0.0 W05-W17 N23-N25
      8  854 1356 (N11E18 L= 16)   342    102        0 ( 827/1309- 887/1399) -33.7 E20-E14 N09-N13
      9 1771 1283 (N32W45 L= 79)   497    105        2 (1732/1242-1801/1305) -15.9 W41-W48 N30-N33
     10  912 1159 (N01E10 L= 24)   219     95        1 ( 854/1144- 959/1178)  31.0 E12-E07 N00-N03
     11 1750 1029 (N16W47 L= 81)   319     48        0 (1725/1012-1778/1044)  11.3 W44-W49 N14-N16
     12  887  488 (S39W07 L= 40)   308     53        1 ( 852/ 467- 917/ 513)   9.5 W03-W09 S40-S37
  

    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