Color  xc   yc  (LatLon LCarr)  Area Length Branches (  x0/y0     x1/y1  )  Tilt Location
      1 1049 1757 (N52W31 L=290)   647    103        2 (1022/1719-1084/1790)  16.7 W27-W37 N48-N53
      2  930 1685 (N51W16 L=275)   327     58        0 ( 898/1662- 956/1702) -63.4 W13-W18 N48-N52
      3 1676 1633 (N22W74 L=333)   817     71        0 (1648/1611-1697/1661) -26.6 W69-W79 N20-N23
      4  553 1361 (N38E25 L=234)  1090    209        3 ( 466/1294- 646/1397) -31.0 E32-E17 N32-N41
      5 1066 1341 (N24W11 L=270)  1175    357        1 ( 915/1313-1245/1368) -18.4 W01-W22 N20-N27
      6 1767 1140 (S06W51 L=310)  1060    151        4 (1734/1071-1794/1220)  45.0 W47-W56 S09-S02
      7  580 1078 (N20E26 L=233)   913    243        0 ( 488/ 992- 667/1140) -50.2 E32-E21 N13-N25
      8  200  967 (N20E63 L=196)  1242    129        2 ( 159/ 905- 238/1013) -41.2 E68-E58 N15-N22
      9 1679  635 (S35W36 L=295)  6574    896       14 (1467/ 381-1859/1006)  35.5 W15-W58 S48-S17
     10 1298  908 (S06W13 L=272)   247     65        0 (1278/ 874-1318/ 947) -90.0 W11-W13 S08-S03
     11 1241  870 (S07W09 L=268)   264    100        0 (1192/ 858-1291/ 881) -26.6 W05-W11 S08-S05
     12 1253  648 (S20W04 L=263)   760    310        1 (1166/ 532-1299/ 771)  47.5 E03-W08 S25-S13
     13 1350  312 (S49W03 L=262)   541     83        2 (1322/ 273-1372/ 349)  39.8 E01-W05 S52-S45
     14 1267  264 (S51E07 L=252)   620     69        1 (1233/ 239-1302/ 293) -50.2 E08-E03 S53-S47
     15  946  244 (S41E32 L=227)   572     91        1 ( 905/ 223- 984/ 270)   6.3 E36-E27 S42-S39
     16 1161  214 (S52E19 L=240)   571     57        0 (1128/ 191-1199/ 234)   7.1 E22-E14 S53-S49
  

    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