Color  xc   yc  (LatLon LCarr)  Area Length Branches (  x0/y0     x1/y1  )  Tilt Location
      1  796 1797 (N34E43 L=299)   554     47        0 ( 774/1765- 818/1831) -56.3 E46-E40 N31-N37
      2  697 1716 (N26E44 L=299)   381     67        1 ( 677/1684- 718/1748) -68.2 E45-E41 N23-N28
      3 1318 1656 (N40E00 L=342)  2064    275        7 (1235/1557-1399/1742) -29.4 E08-W08 N35-N45
      4  989 1536 (N21E16 L=326)   633    211        0 ( 906/1475-1071/1561)   0.0 E21-E09 N18-N23
      5  763 1419 (N09E26 L=317)   416    124        1 ( 744/1361- 780/1473) -50.2 E27-E22 N05-N11
      6  584 1378 (N02E35 L=307)   339     87        0 ( 567/1340- 611/1418) -45.0 E37-E32 N00-N04
      7 1696 1333 (N31W36 L= 18)  1406    217        3 (1624/1262-1762/1402) -14.9 W28-W43 N27-N33
      8 1408 1166 (N11W18 L=  0)   724    178        0 (1385/1084-1446/1247) -84.3 W16-W20 N06-N16
      9  462 1176 (S12E38 L=304)   532    213        2 ( 383/1153- 569/1200)  21.0 E43-E30 S14-S09
     10 1601 1090 (N13W32 L= 15)   848    236        1 (1545/ 991-1656/1192) -30.3 W26-W38 N08-N16
     11  194 1105 (S21E63 L=280)   659    114        0 ( 150/1059- 223/1152)  33.7 E68-E59 S23-S17
     12  633  976 (S19E22 L=320)   349     74        0 ( 590/ 956- 679/ 999)  38.7 E24-E19 S21-S17
     13  440  918 (S27E35 L=307)   231     79        0 ( 404/ 901- 472/ 942)  31.0 E37-E32 S28-S25
     14  576  735 (S35E22 L=321)   732    161        2 ( 506/ 704- 649/ 772)  35.0 E27-E16 S38-S31
     15 1733  542 (S11W66 L= 48)   772    110        0 (1685/ 508-1778/ 575)  54.5 W62-W69 S14-S07
     16 1341  494 (S27W37 L= 20)  1264    212        2 (1254/ 451-1437/ 559)  63.4 W32-W40 S32-S20
     17  588  497 (S51E17 L=326)   257     42        0 ( 569/ 476- 611/ 520)  51.3 E18-E14 S53-S48
     18 1514  326 (S30W64 L= 47)   702    124        1 (1460/ 298-1578/ 354)  36.9 W60-W68 S32-S26
     19 1117  337 (S44W34 L= 16)   428     72        0 (1080/ 322-1158/ 352)  21.8 W30-W36 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