Color  xc   yc  (LatLon LCarr)  Area Length Branches (  x0/y0     x1/y1  )  Tilt Location
      1 1135 1783 (N50E15 L= 50)  4210    372        6 (1051/1626-1218/1913) -35.7 E37-E03 N38-N63
      2  379 1536 (N15E59 L=  6)  2024    225        6 ( 334/1447- 418/1649) -45.0 E69-E53 N08-N21
      3  602 1224 (N01E29 L= 36)  1289    275        2 ( 532/1119- 685/1325)  69.0 E32-E26 S06-N07
      4 1365 1264 (N19W15 L= 80)   204     69        0 (1332/1248-1394/1289)   0.0 W12-W17 N18-N19
      5  801  995 (S09E12 L= 53)   443    144        0 ( 743/ 950- 864/1035)  49.4 E15-E09 S12-S05
      6 1772  792 (N01W56 L=121)  3076    263        6 (1707/ 684-1831/ 930) -32.9 W47-W66 S04-N08
      7 1429  865 (S03W27 L= 92)   570     98        1 (1389/ 821-1470/ 908)  68.2 W25-W29 S05-N00
      8  430  684 (S36E35 L= 30)   460    156        1 ( 380/ 624- 454/ 716) -39.8 E39-E32 S39-S33
      9 1742  646 (S08W59 L=124)   526     54        1 (1713/ 623-1768/ 673)  -7.1 W55-W63 S08-S06
     10 1173  172 (S54W51 L=116)   462     69        0 (1150/ 145-1194/ 201) -33.7 W46-W57 S56-S50
  

    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