Color  xc   yc  (LatLon LCarr)  Area Length Branches (  x0/y0     x1/y1  )  Tilt Location
      1 1231 1572 (N36E01 L=257)   263     52        1 (1209/1550-1251/1601) -71.6 E02-W00 N34-N37
      2  778 1511 (N20E28 L=230)   521    157        1 ( 763/1441- 796/1580) -66.0 E30-E25 N15-N24
      3  675 1518 (N18E35 L=223)   242     81        0 ( 664/1480- 689/1552) -53.1 E36-E32 N15-N19
      4 1087 1471 (N25E07 L=250)   426     93        0 (1048/1428-1128/1510) -23.2 E10-E03 N22-N26
      5  572 1463 (N12E39 L=218)   264     54        1 ( 553/1440- 587/1493) -45.0 E41-E37 N11-N14
      6  647 1368 (N08E31 L=227)   351     67        2 ( 633/1329- 665/1407) -68.2 E32-E29 N05-N10
      7  220 1177 (S11E60 L=198)  1385    168        1 ( 180/1104- 258/1260)  68.2 E63-E58 S16-S06
      8 1334 1181 (N13W14 L=272)   247     86        0 (1293/1162-1374/1196)  36.9 W12-W16 N11-N14
      9  572 1030 (S13E27 L=230)   563    137        3 ( 514/ 997- 628/1063)  23.2 E31-E24 S15-S09
     10  184  995 (S23E60 L=198)   773    110        0 ( 156/ 944- 201/1032) -56.3 E63-E57 S26-S20
     11 1205  753 (S14W17 L=275)   387     82        1 (1169/ 732-1241/ 783)  53.1 W14-W19 S15-S11
     12 1405  711 (S12W31 L=288)   243     76        0 (1378/ 681-1426/ 749)  76.0 W29-W31 S13-S09
     13  555  700 (S33E23 L=234)   631    116        0 ( 529/ 646- 595/ 737)  74.1 E25-E20 S37-S30
     14  343  693 (S39E42 L=216)   478     86        0 ( 316/ 664- 378/ 718)  33.7 E45-E38 S41-S36
     15  351  591 (S46E43 L=215)   492     88        0 ( 323/ 553- 379/ 640) -51.3 E45-E40 S48-S43
     16 1109  540 (S29W19 L=276)   355     83        1 (1060/ 523-1154/ 561)  18.4 W15-W21 S31-S27
     17 1631  470 (S19W59 L=316)   519     53        2 (1602/ 449-1659/ 490)  36.9 W56-W61 S20-S17
     18  877  443 (S42W07 L=264)   277     41        1 ( 851/ 426- 905/ 465)  33.7 W04-W08 S43-S40
  

    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