Color  xc   yc  (LatLon LCarr)  Area Length Branches (  x0/y0     x1/y1  )  Tilt Location
      1 1168 1582 (N40W07 L=340)   502    105        1 (1118/1554-1219/1598)   7.1 W03-W11 N37-N41
      2  671 1481 (N29E29 L=303)   517    136        0 ( 616/1459- 741/1513) -10.3 E34-E23 N27-N30
      3 1475 1424 (N30W31 L=  3)   332     84        1 (1437/1406-1514/1447)   0.0 W27-W33 N29-N31
      4  527 1299 (N16E36 L=297)   964    287        1 ( 453/1233- 629/1400) -27.9 E44-E27 N12-N21
      5 1003 1303 (N19E03 L=329)   585    199        0 ( 919/1242-1078/1344) -26.6 E08-W02 N16-N21
      6  555  761 (S17E30 L=303)   237     85        0 ( 537/ 720- 572/ 804)  59.0 E31-E28 S19-S14
      7 1481  695 (S16W33 L=  5)   370     62        1 (1455/ 666-1504/ 730)  90.0 W30-W34 S17-S13
      8 1715  646 (S18W55 L= 27)   351     44        1 (1690/ 621-1734/ 665)  63.4 W52-W56 S19-S16
      9  308  540 (S35E61 L=271)   845    112        1 ( 259/ 517- 355/ 571)  12.3 E75-E52 S37-S32
  

    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