Color  xc   yc  (LatLon LCarr)  Area Length Branches (  x0/y0     x1/y1  )  Tilt Location
      1  509 1744 (N24E60 L=193)  1515    188        2 ( 453/1659- 552/1807) -47.5 E69-E56 N17-N29
      2 1676 1484 (N39W33 L=287)   633     74        2 (1641/1449-1703/1519) -53.1 W29-W36 N37-N42
      3 1818 1341 (N35W49 L=303)  1096    179        0 (1734/1287-1883/1385)   0.0 W39-W58 N33-N36
      4  213 1281 (S09E63 L=191)  3980    359       14 ( 115/1134- 294/1374) -28.4 E81-E53 S17-S03
      5  171 1233 (S12E67 L=187)  1117    134        1 ( 143/1180- 204/1292)  54.5 E70-E64 S15-S08
      6  494 1170 (S11E35 L=219)   284     96        0 ( 472/1134- 523/1209) -38.7 E37-E32 S12-S08
      7 1385  682 (S16W30 L=284)   213    103        0 (1360/ 633-1417/ 723) -80.5 W29-W31 S19-S13
      8 1525  609 (S16W42 L=296)   359     76        1 (1498/ 582-1545/ 652) -31.0 W40-W45 S17-S13
      9  850  433 (S47W07 L=261)  1798    403        3 ( 684/ 359-1040/ 503)   5.4 E08-W24 S49-S43
  

    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