Color  xc   yc  (LatLon LCarr)  Area Length Branches (  x0/y0     x1/y1  )  Tilt Location
      1  797 1653 (N52E06 L=248)   502     67        0 ( 768/1618- 820/1681) -63.4 E08-E03 N48-N54
      2 1815 1319 (N08W65 L=319)   752     73        1 (1776/1291-1854/1346) -18.4 W60-W70 N05-N09
      3  250 1081 (N19E56 L=197)  1242    194        2 ( 184/1049- 334/1129)   6.7 E65-E47 N17-N21
      4 1268  920 (S04W13 L=267)   428    119        1 (1216/ 898-1336/ 937) -14.0 W09-W17 S04-S02
      5  322  841 (N04E51 L=203)   372     41        0 ( 297/ 815- 350/ 865)  31.0 E53-E48 N02-N05
      6 1383  534 (S31W16 L=270)  2959    578        4 (1160/ 396-1564/ 718)  22.8 E01-W30 S37-S21
      7  312  614 (S10E59 L=194)   476     50        1 ( 289/ 594- 338/ 637)   9.5 E62-E56 S11-S08
      8  582  535 (S18E38 L=215)   525     82        0 ( 561/ 492- 610/ 574)  45.0 E40-E34 S20-S15
      9  822  373 (S32E27 L=226)   458     88        1 ( 793/ 342- 858/ 415) -81.9 E28-E24 S35-S28
     10  577  343 (S31E49 L=204)  1081    176        2 ( 551/ 272- 624/ 414)  34.5 E58-E40 S37-S26
  

    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