Color  xc   yc  (LatLon LCarr)  Area Length Branches (  x0/y0     x1/y1  )  Tilt Location
      1  379 1542 (N16E55 L= 34)   578    143        1 ( 362/1501- 398/1591) -50.2 E58-E52 N12-N18
      2 1718 1528 (N40W43 L=133)   428     95        0 (1697/1500-1742/1551)  50.2 W40-W46 N37-N43
      3 1542 1422 (N27W27 L=116)   967    214        1 (1436/1366-1626/1483)  -7.6 W18-W33 N24-N29
      4 1483 1142 (N09W25 L=115)   847    203        3 (1448/1023-1511/1234) -72.9 W22-W28 N01-N14
      5 1701 1069 (N09W42 L=131)   522    144        2 (1667/1011-1739/1117) -63.4 W38-W45 N05-N11
      6  267 1058 (S15E52 L= 37)   379     73        2 ( 238/1034- 297/1083)  -9.5 E54-E48 S16-S13
      7 1860  889 (N03W61 L=150)   706    174        0 (1824/ 836-1911/ 922)  45.0 W57-W67 N00-N06
      8 1103  771 (S20W09 L= 99)   702    197        0 (1012/ 714-1203/ 851)  50.7 W04-W13 S24-S13
      9 1695  512 (S22W59 L=149)   537     84        1 (1671/ 477-1724/ 543) -76.0 W56-W61 S23-S19
     10 1070  367 (S46W20 L=109)   592     99        1 (1025/ 340-1125/ 394)   0.0 W15-W25 S47-S44
     11  863  204 (S67W11 L=101)   525     77        1 ( 838/ 177- 883/ 226) -21.8 W06-W16 S69-S63
  

    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