Color  xc   yc  (LatLon LCarr)  Area Length Branches (  x0/y0     x1/y1  )  Tilt Location
      1 1082 1653 (N41W26 L=257)   543    102        2 (1040/1607-1123/1690)  20.6 W22-W30 N38-N42
      2 1585 1604 (N21W58 L=289)  1283    190        3 (1490/1569-1654/1648)  -5.4 W48-W69 N18-N23
      3 1539 1427 (N13W43 L=274)  1203    210        1 (1458/1350-1601/1492) -76.0 W40-W45 N07-N19
      4  817 1345 (N30E03 L=227)   747    196        2 ( 732/1297- 907/1376) -45.0 E08-W01 N24-N33
      5  733 1227 (N25E12 L=219)   241     99        0 ( 706/1180- 772/1258) -71.6 E13-E10 N21-N27
      6 1262 1122 (N06W16 L=246)  1215    351        6 (1200/ 956-1343/1246)  32.0 W07-W23 S01-N11
      7 1421 1191 (N05W27 L=258)   394     95        1 (1381/1154-1463/1235) -66.8 W25-W28 N01-N08
      8 1927 1041 (S21W67 L=298)   894     82        2 (1891/1017-1955/1072)   0.0 W60-W74 S22-S19
      9 1657  950 (S15W36 L=266)   655    147        4 (1583/ 912-1708/ 976) -16.7 W30-W40 S19-S12
     10  534  873 (N10E33 L=198)   279     40        1 ( 511/ 853- 555/ 895)  14.0 E34-E30 N08-N11
     11 1868  838 (S30W56 L=286)   887     73        3 (1831/ 814-1908/ 863) -14.0 W50-W62 S32-S28
     12  783  735 (S03E21 L=210)   253     90        1 ( 763/ 692- 804/ 782)  78.7 E21-E20 S05-N00
     13  299  715 (N06E56 L=175)  1044    215        3 ( 225/ 682- 402/ 740)  -6.3 E65-E47 N02-N08
     14  946  630 (S13E15 L=216)   206     53        2 ( 923/ 606- 964/ 653) -71.6 E16-E13 S14-S11
     15  482  465 (S12E52 L=179)  9624      0       32 ( 316/ 303- 640/ 615) -65.0 E79-E36 S23-N00
     16  764  393 (S23E35 L=195)  1393    327        3 ( 708/ 273- 836/ 563)  37.7 E46-E24 S30-S13
     17  652  515 (S13E37 L=194)   238     44        2 ( 630/ 496- 668/ 539)  33.7 E38-E35 S14-S11
     18 1575  446 (S47W22 L=253)   410     60        0 (1562/ 411-1593/ 479) -90.0 W20-W24 S49-S43
     19  895  269 (S35E35 L=196)  1079    150        3 ( 825/ 225- 958/ 318)   8.1 E41-E27 S37-S32
  

    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