Color  xc   yc  (LatLon LCarr)  Area Length Branches (  x0/y0     x1/y1  )  Tilt Location
      1  362 1565 (N29E64 L=212)  2005    293        4 ( 325/1461- 401/1665) -24.1 E86-E53 N22-N35
      2  571 1503 (N27E39 L=238)   302     35        1 ( 549/1482- 599/1524)  56.3 E40-E36 N25-N28
      3  896 1236 (N13E10 L=266)   239    128        0 ( 843/1214- 955/1251)   0.0 E13-E06 N12-N14
      4 1540  973 (N03W34 L=310)   625    160        1 (1497/ 911-1586/1031)  54.5 W31-W36 S01-N06
      5 1028  983 (S01W01 L=277)   462    124        1 ( 973/ 960-1063/1028)  14.0 E02-W02 S02-N01
      6  170  854 (S18E68 L=208)  2156    166        4 ( 145/ 771- 197/ 935)  70.0 E73-E63 S23-S12
      7  906  828 (S12E06 L=271)   330     70        0 ( 865/ 807- 941/ 845) -11.3 E08-E03 S12-S10
      8 1716  680 (S14W54 L=331)  2060    250        5 (1671/ 598-1767/ 781) -76.0 W50-W61 S19-S07
      9 1089  718 (S17W07 L=283)   510    165        0 (1014/ 696-1152/ 763) -11.3 W01-W11 S18-S14
     10 1486  671 (S16W35 L=311)   711    225        0 (1427/ 631-1546/ 738)  49.4 W29-W38 S19-S11
     11  318  709 (S25E51 L=225)   574    111        0 ( 267/ 690- 363/ 735)  14.0 E57-E45 S26-S23
     12  915  505 (S33E02 L=274)   544    138        1 ( 844/ 483- 970/ 529)  14.0 E07-W01 S34-S30
  

    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