Color  xc   yc  (LatLon LCarr)  Area Length Branches (  x0/y0     x1/y1  )  Tilt Location
      1 1561 1613 (N21W55 L=289)   427     70        0 (1519/1601-1601/1626) -26.6 W51-W58 N19-N23
      2  456 1392 (N41E30 L=205)   366    115        0 ( 413/1356- 507/1418) -36.9 E33-E25 N37-N43
      3 1373 1191 (N05W24 L=258)  1328    414        2 (1222/1070-1536/1279)   4.8 W12-W36 N01-N07
      4 1022 1183 (N14W04 L=238)   239     56        2 ( 994/1156-1051/1208)   0.0 W02-W05 N12-N14
      5  496 1125 (N23E30 L=205)   241     97        2 ( 458/1097- 527/1157)  21.8 E32-E27 N22-N24
      6 1926  995 (S23W66 L=300)  1290    125        2 (1893/ 948-1958/1050) -63.4 W59-W72 S25-S19
      7  750  983 (N09E16 L=218)   754    160        2 ( 700/ 927- 796/1047)  23.2 E19-E12 N05-N12
      8 1425  948 (S10W20 L=254)   337    130        0 (1393/ 900-1473/1007) -76.0 W19-W22 S13-S05
      9 1808  820 (S29W48 L=282)  1082    194        1 (1754/ 762-1864/ 900)  21.0 W42-W55 S31-S26
     10  474  598 (S06E45 L=189)  1644    342        4 ( 409/ 476- 545/ 749)  82.9 E50-E42 S13-N03
     11  835  647 (S11E21 L=214)   243     43        1 ( 818/ 615- 853/ 671)  76.0 E21-E19 S13-S09
     12 1617  602 (S37W28 L=262)  1302    189        3 (1536/ 541-1685/ 661)   4.1 W20-W34 S39-S33
     13  382  491 (S10E59 L=175)  1865    328        3 ( 349/ 359- 430/ 626)  32.6 E75-E50 S18-S02
     14 1092  535 (S25E10 L=225)   734    186        0 (1053/ 458-1116/ 620)  48.4 E14-E06 S28-S19
     15  600  532 (S12E39 L=195)   454    142        0 ( 566/ 482- 653/ 562)   0.0 E42-E35 S15-S10
     16  621  422 (S19E43 L=191)   815    164        2 ( 569/ 364- 684/ 479)  83.7 E45-E40 S24-S15
     17  531  376 (S20E53 L=181)  1661    262        3 ( 454/ 312- 637/ 446) -42.0 E59-E45 S26-S14
     18  767  278 (S33E43 L=191)   493    126        1 ( 697/ 253- 811/ 305) -12.5 E47-E37 S34-S30
  

    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