Color  xc   yc  (LatLon LCarr)  Area Length Branches (  x0/y0     x1/y1  )  Tilt Location
      1 1515 1730 (N35W59 L=165)   882    175        1 (1449/1708-1566/1744) -21.8 W53-W63 N32-N37
      2 1019 1578 (N41W15 L=121)   814    188        5 ( 972/1520-1080/1624) -76.0 W11-W18 N35-N44
      3 1468 1288 (N12W33 L=139)  2914    767        8 (1158/1079-1688/1468) -54.8 W19-W44 S05-N29
      4  447 1283 (N34E35 L= 71)   322     99        0 ( 417/1252- 487/1314) -33.7 E37-E31 N31-N36
      5 1635  861 (S17W35 L=141)   805    135        2 (1573/ 819-1676/ 903) -40.6 W30-W38 S20-S14
      6  525  513 (S14E45 L= 61)  1348    305        2 ( 497/ 397- 561/ 615)  49.4 E51-E38 S22-S08
      7 1350  449 (S37W08 L=114)  1330    277        6 (1301/ 388-1404/ 519)  73.3 W05-W12 S42-S31
      8  828  397 (S27E28 L= 78)   289     74        1 ( 788/ 379- 863/ 419) -36.9 E30-E25 S28-S25
      9  658  390 (S25E40 L= 66)   294     76        2 ( 626/ 372- 687/ 406)   0.0 E43-E38 S25-S23
  

    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