Color  xc   yc  (LatLon LCarr)  Area Length Branches (  x0/y0     x1/y1  )  Tilt Location
      1  477 1733 (N30E67 L=345)   987    113        0 ( 457/1675- 492/1785) -21.8 E82-E61 N25-N33
      2 1602 1505 (N41W33 L= 85)   444     81        0 (1563/1473-1635/1529) -15.9 W28-W35 N39-N42
      3  798 1168 (N01E16 L= 36)  1175    295        0 ( 721/1054- 869/1285)  75.1 E18-E13 S06-N09
      4 1149 1253 (N14W02 L= 55)   233     54        0 (1132/1225-1176/1273)  76.0 W01-W03 N11-N15
      5  196 1237 (S06E65 L=348)   478     60        0 ( 176/1214- 213/1267)  76.0 E67-E61 S07-S03
      6 1519 1205 (N19W27 L= 79)   367     99        0 (1467/1172-1561/1234)   0.0 W22-W30 N17-N19
      7 1010  922 (S08W01 L= 54)   208     71        0 ( 985/ 891-1042/ 942)  63.4 W00-W02 S10-S06
      8 1595  793 (S03W41 L= 93)   799    130        1 (1540/ 752-1631/ 851)  56.3 W37-W42 S06-N00
      9 1872  790 (N03W68 L=121)  1012    105        1 (1828/ 752-1905/ 821)  -3.8 W61-W76 N01-N04
  

    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