Color  xc   yc  (LatLon LCarr)  Area Length Branches (  x0/y0     x1/y1  )  Tilt Location
      1 1691 1430 (N26W55 L=346)   951    154        2 (1642/1366-1736/1499) -71.6 W52-W58 N21-N30
      2 1113 1426 (N29W08 L=299)   648    130        4 (1081/1379-1145/1491)  60.3 W05-W11 N26-N33
      3  294 1382 (N28E59 L=231)  1476    228        2 ( 251/1303- 340/1455) -84.3 E66-E55 N22-N32
      4 1193 1329 (N23W13 L=304)   496    121        0 (1134/1310-1248/1351) -20.6 W08-W16 N21-N24
      5  957  867 (S05E05 L=286)   216     50        1 ( 937/ 839- 982/ 888)  33.7 E06-E03 S06-S03
      6 1749  853 (S11W51 L=342)   587     58        2 (1719/ 823-1774/ 881) -63.4 W48-W53 S12-S08
      7 1068  804 (S09W02 L=293)   355     68        2 (1042/ 771-1095/ 833) -53.1 W00-W03 S11-S07
      8  886  780 (S10E10 L=281)   259     78        1 ( 869/ 748- 908/ 824) -78.7 E10-E08 S12-S07
      9 1339  744 (S14W19 L=310)   269     77        0 (1324/ 705-1350/ 786)  78.7 W18-W20 S16-S11
     10 1697  607 (S27W51 L=342)   843    132        0 (1668/ 579-1742/ 649)   8.1 W47-W55 S28-S23
     11 1693  496 (S35W57 L=348)   480     40        0 (1664/ 482-1719/ 512)  -7.1 W52-W60 S36-S34
  

    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