Color  xc   yc  (LatLon LCarr)  Area Length Branches (  x0/y0     x1/y1  )  Tilt Location
      1 1322 1731 (N49W40 L= 48)   482     67        2 (1292/1698-1346/1759)  78.7 W38-W43 N45-N50
      2 1269 1620 (N41W28 L= 35)   482     91        1 (1229/1582-1301/1643) -45.0 W24-W30 N37-N42
      3 1599 1605 (N34W58 L= 65)   428     55        1 (1577/1579-1620/1638)  53.1 W55-W60 N32-N36
      4 1708 1490 (N25W62 L= 69)   766    108        1 (1664/1451-1746/1521) -45.0 W57-W64 N22-N27
      5  354 1425 (N35E52 L=315)   450     73        1 ( 319/1408- 388/1447) -11.3 E57-E47 N34-N36
      6 1057 1417 (N29W06 L= 14)   297    121        1 (1000/1391-1108/1443) -29.7 W02-W09 N27-N31
      7 1580 1058 (N01W36 L= 44)   228     51        1 (1558/1032-1598/1086)  33.7 W34-W37 N00-N02
      8  516  829 (S03E35 L=333)   232     76        0 ( 485/ 799- 540/ 857) -53.1 E36-E33 S04-N00
      9 1733  794 (S17W48 L= 56)   467    130        1 (1703/ 750-1752/ 830)  90.0 W45-W50 S20-S14
     10  391  720 (S09E47 L=321)   820    188        3 ( 355/ 641- 422/ 792) -84.3 E49-E43 S14-S04
     11 1013  702 (S14E04 L=  4)   279    180        0 ( 966/ 677-1070/ 754) -45.0 E06-W00 S16-S10
     12 1117  640 (S19W02 L= 10)   522    179        1 (1041/ 602-1157/ 708)  39.8 E02-W04 S21-S15
     13 1669  615 (S29W44 L= 52)   825    220        1 (1576/ 582-1763/ 651)  -3.0 W35-W54 S30-S27
     14 1261  532 (S29W11 L= 19)  1387    423        1 (1122/ 433-1442/ 649)  24.6 W00-W24 S34-S22
     15  832  548 (S23E18 L=350)   406    136        0 ( 801/ 493- 884/ 601)  45.0 E20-E13 S26-S19
     16 1281  425 (S37W13 L= 20)   538    158        3 (1212/ 405-1366/ 452)   0.0 W07-W19 S38-S35
     17  938  293 (S44E17 L=351)  1436    210        2 ( 887/ 217-1015/ 349)  36.9 E24-E08 S51-S39
     18 1030  215 (S54E12 L=355)   364     64        0 (1005/ 197-1063/ 236) -51.3 E13-E09 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