Color  xc   yc  (LatLon LCarr)  Area Length Branches (  x0/y0     x1/y1  )  Tilt Location
      1  487 1497 (N31E44 L= 37)   601     87        0 ( 437/1478- 539/1512)   0.0 E49-E39 N29-N32
      2  905 1363 (N24E09 L= 72)   529    198        0 ( 827/1348-1022/1378)   0.0 E14-E01 N23-N24
      3 1461 1199 (N15W28 L=109)   309    122        0 (1400/1181-1502/1218)   0.0 W24-W31 N13-N15
      4  146  872 (S11E73 L=  8)  2404    208        5 ( 109/ 792- 173/ 950) -48.4 E80-E68 S15-S06
      5 1390  844 (S07W24 L=105)   259     68        0 (1361/ 812-1421/ 868)  45.0 W21-W25 S09-S05
      6 1631  750 (S13W43 L=124)   669    179        0 (1573/ 687-1699/ 816)  45.0 W39-W47 S17-S09
      7  351  723 (S19E48 L= 32)  1099    217        1 ( 322/ 624- 404/ 809) -80.5 E50-E45 S25-S13
      8 1635  695 (S17W44 L=125)   326     58        0 (1614/ 671-1650/ 726)  76.0 W42-W45 S18-S14
      9 1039  502 (S31W03 L= 83)   573    129        1 (1020/ 438-1063/ 562) -71.6 W01-W04 S35-S26
     10 1092  374 (S41W08 L= 89)   332     34        0 (1070/ 348-1112/ 401)  63.4 W06-W09 S43-S38
  

    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