Color  xc   yc  (LatLon LCarr)  Area Length Branches (  x0/y0     x1/y1  )  Tilt Location
      1  557 1545 (N12E43 L=189)  1076    299        0 ( 510/1409- 630/1665) -80.5 E45-E40 N03-N21
      2 1317 1231 (N12W11 L=244)  2015    248        6 (1216/1047-1451/1454) -39.3 W00-W22 N04-N22
      3 1473 1424 (N28W17 L=249)   228     38        0 (1450/1400-1491/1447) -18.4 W15-W18 N26-N29
      4  663 1332 (N02E28 L=204)  1466    407        5 ( 480/1249- 797/1373)   0.0 E41-E18 S01-N04
      5  780 1155 (S06E17 L=215)   510    144        1 ( 759/1083- 800/1234) -58.0 E19-E14 S09-S01
      6 1784 1002 (N13W49 L=281)  2865    426        2 (1747/ 805-1819/1216) -57.5 W42-W58 N03-N25
      7  598 1086 (S13E27 L=206)   925    237        1 ( 485/1045- 706/1140)  34.7 E33-E20 S17-S08
      8  990  726 (S25W06 L=238)  1118    372        0 ( 807/ 696-1168/ 766)  30.1 E04-W16 S30-S18
      9  456  714 (S39E31 L=201)  1578    377        4 ( 312/ 690- 602/ 741)  15.1 E45-E19 S43-S35
  

    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