Color  xc   yc  (LatLon LCarr)  Area Length Branches (  x0/y0     x1/y1  )  Tilt Location
      1 1273 1718 (N44W40 L= 84)   573    131        1 (1235/1682-1327/1748)   5.2 W34-W45 N42-N45
      2 1107 1357 (N25W12 L= 56)   301     98        1 (1073/1323-1150/1382)  18.4 W09-W15 N23-N26
      3  254 1342 (N39E55 L=349)   441     70        1 ( 232/1312- 277/1370)  36.9 E56-E52 N36-N40
      4  722 1339 (N31E14 L= 30)   222    105        0 ( 668/1330- 769/1347) -15.9 E17-E10 N30-N32
      5  294 1087 (N22E50 L=354)   750    185        3 ( 232/1061- 350/1125) -10.3 E55-E44 N20-N24
      6  142  933 (N14E70 L=334)   868    145        1 ( 114/ 892- 163/ 974) -45.0 E72-E66 N11-N16
      7  532  951 (N11E32 L= 12)   270     68        0 ( 498/ 934- 563/ 965) -21.8 E34-E29 N09-N12
      8 1335  347 (S45W07 L= 51)   384     71        0 (1306/ 331-1373/ 367) -45.0 W05-W10 S47-S42
  

    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