Color  xc   yc  (LatLon LCarr)  Area Length Branches (  x0/y0     x1/y1  )  Tilt Location
      1 1827 1283 (N16W64 L=119)  1138    185        1 (1795/1211-1856/1351)  66.0 W58-W68 N11-N20
      2  502 1277 (N20E35 L= 20)   267     18        1 ( 486/1254- 522/1299)  90.0 E36-E33 N18-N21
      3 1133 1042 (N05W07 L= 62)   196     59        0 (1102/1026-1175/1055)  11.3 W04-W09 N03-N05
      4  307  979 (N02E50 L=  6)   409     97        0 ( 261/ 957- 350/1008) -18.4 E54-E45 N00-N03
      5 1151  807 (S10W07 L= 63)  1385    571        1 ( 926/ 655-1268/ 987)  45.0 E07-W14 S18-N01
      6 1734  935 (S05W49 L=104)   541     76        0 (1695/ 912-1767/ 970)  26.6 W45-W51 S05-S02
      7  785  914 (S02E15 L= 41)   376     86        1 ( 753/ 875- 823/ 948)  53.1 E17-E12 S04-N00
      8  467  855 (S06E37 L= 19)   518    148        0 ( 409/ 820- 504/ 911)  29.7 E41-E34 S07-S02
      9  255  786 (S10E57 L=358)   925    140        2 ( 224/ 722- 291/ 848)  58.0 E60-E52 S14-S06
     10  200  734 (S14E66 L=349)  3375    205        7 ( 156/ 630- 248/ 828) -80.5 E74-E61 S20-S08
     11  290  676 (S17E56 L=359)   402     67        0 ( 266/ 646- 316/ 707)  33.7 E60-E52 S19-S15
  

    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