Color  xc   yc  (LatLon LCarr)  Area Length Branches (  x0/y0     x1/y1  )  Tilt Location
      1  961 1796 (N58W29 L=223)   382     62        2 ( 937/1763- 980/1827)  20.0 W23-W34 N55-N59
      2  620 1627 (N56E15 L=179)   298     65        1 ( 596/1610- 652/1653)  -8.1 E17-E10 N54-N58
      3  340 1242 (N35E45 L=149)  1927    432        4 ( 201/1153- 517/1306)  -3.4 E62-E28 N30-N38
      4  820 1286 (N27E06 L=188)   321    101        0 ( 766/1271- 871/1300) -15.9 E09-E02 N25-N28
      5 1639 1184 (S01W42 L=235)   394     91        1 (1588/1166-1683/1196) -15.9 W38-W45 S02-N01
      6 1234 1084 (N05W13 L=207)   374    158        0 (1198/1022-1288/1146)  29.1 W09-W18 N02-N07
      7 1735  880 (S21W43 L=236)  1279    336        1 (1611/ 766-1839/ 984)  10.3 W30-W54 S24-S17
      8 1262  918 (S05W11 L=204)   356     75        0 (1239/ 879-1289/ 947) -90.0 W09-W12 S07-S02
      9 1791  791 (S29W48 L=241)   886    193        0 (1754/ 722-1817/ 889)  59.0 W43-W50 S32-S22
     10 1178  613 (S21E01 L=192)  1207    446        2 (1018/ 495-1290/ 766)  22.2 E14-W08 S24-S13
     11 1502  663 (S27W20 L=214)   457    109        0 (1461/ 631-1544/ 707)  15.9 W17-W24 S28-S24
     12 1441  425 (S42W12 L=206)  1651    396        0 (1294/ 319-1560/ 562)  21.8 E02-W23 S46-S35
  

    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