Color  xc   yc  (LatLon LCarr)  Area Length Branches (  x0/y0     x1/y1  )  Tilt Location
      1  731 1866 (N38E58 L=136)   705     57        1 ( 700/1844- 765/1883)   0.0 E64-E53 N36-N39
      2  809 1626 (N23E32 L=162)   709    207        2 ( 768/1532- 858/1709) -34.7 E37-E24 N18-N27
      3  403 1524 (N07E55 L=138)  2076    338        5 ( 333/1397- 509/1632)  86.4 E61-E51 S01-N15
      4 1108 1536 (N25E10 L=183)   337    118        0 (1092/1476-1128/1588) -56.3 E12-E07 N22-N29
      5 1197 1277 (N13W03 L=196)   532    199        1 (1163/1193-1216/1353) -66.0 W00-W04 N07-N16
      6  408 1067 (S19E39 L=154)  1700    315        2 ( 286/ 949- 505/1170)  50.5 E48-E33 S28-S11
      7  134 1102 (S22E70 L=123)   715     86        2 ( 105/1072- 158/1126) -23.2 E76-E65 S23-S20
      8  507  839 (S30E27 L=166)  2922    671        3 ( 221/ 724- 813/ 927)  28.4 E58-E08 S44-S17
      9 1670  843 (N03W44 L=237)   602    121        1 (1632/ 812-1711/ 887) -15.9 W39-W47 N01-N05
     10 1415  796 (S07W28 L=222)   259     82        0 (1388/ 760-1449/ 826)  78.7 W27-W29 S09-S04
     11 1161  623 (S24W20 L=214)   429    141        0 (1092/ 602-1221/ 650)  45.0 W16-W22 S26-S20
     12 1681  576 (S10W57 L=250)   463     64        1 (1666/ 529-1697/ 608) -51.3 W54-W59 S12-S07
     13 1618  526 (S15W55 L=248)   498     84        1 (1593/ 491-1645/ 565)  90.0 W52-W56 S16-S11
     14  560  477 (S52E17 L=176)   373    115        0 ( 529/ 440- 607/ 521) -18.4 E21-E11 S54-S49
     15 1642  441 (S18W63 L=257)   602     60        1 (1619/ 407-1674/ 464)  90.0 W59-W66 S19-S15
     16 1017  296 (S49W30 L=224)  1928    345        4 ( 860/ 225-1137/ 387)   1.8 W12-W45 S50-S45
  

    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