Color  xc   yc  (LatLon LCarr)  Area Length Branches (  x0/y0     x1/y1  )  Tilt Location
      1  493 1731 (N56E46 L= 80)   783     61        0 ( 471/1701- 524/1763) -21.3 E56-E38 N52-N59
      2 1344 1609 (N32W31 L=158)   328     84        1 (1302/1575-1382/1631) -31.0 W28-W33 N29-N34
      3  229 1403 (N33E62 L= 64)   869     97        1 ( 199/1350- 257/1446) -31.6 E70-E56 N28-N36
      4  764 1351 (N22E13 L=113)   541    163        1 ( 688/1316- 825/1392)  24.0 E17-E08 N20-N24
      5  563 1299 (N21E27 L= 99)   628    214        2 ( 459/1272- 672/1317) -14.9 E34-E19 N18-N23
      6 1858  944 (S14W60 L=186)  1791    289        0 (1836/ 803-1882/1085)  86.6 W57-W63 S23-S05
      7  301  719 (S10E55 L= 72)   875    170        0 ( 262/ 654- 321/ 799)  45.0 E61-E51 S13-S05
      8  394  576 (S20E51 L= 75)   426     76        0 ( 359/ 546- 432/ 598)  12.5 E56-E47 S21-S18
      9  532  549 (S23E41 L= 85)   315    114        0 ( 492/ 522- 575/ 585) -51.3 E42-E37 S25-S20
     10  337  483 (S25E64 L= 62)   541     55        0 ( 305/ 469- 371/ 498)   0.0 E69-E58 S25-S24
     11  727  347 (S40E35 L= 91)   375     59        1 ( 701/ 327- 768/ 366)   8.1 E37-E30 S41-S38
     12  611  272 (S44E53 L= 74)   480     73        1 ( 578/ 249- 643/ 291)   9.5 E59-E47 S45-S42
  

    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