Color  xc   yc  (LatLon LCarr)  Area Length Branches (  x0/y0     x1/y1  )  Tilt Location
      1  624 1704 (N25E48 L= 75)   413     74        1 ( 586/1688- 660/1720)  36.9 E50-E46 N23-N26
      2  811 1436 (N13E24 L= 99)   736    198        2 ( 791/1348- 844/1531) -51.3 E27-E19 N08-N18
      3  260 1189 (S14E56 L= 67)   855    147        1 ( 229/1112- 292/1266) -84.3 E58-E53 S18-S08
      4  690 1199 (S03E24 L= 99)   309     64        1 ( 667/1171- 718/1234)  90.0 E24-E22 S05-S01
      5 1321 1102 (N07W14 L=137)   259     70        0 (1283/1078-1360/1121)  45.0 W12-W16 N05-N09
      6 1663 1023 (N13W39 L=162)   724    205        2 (1603/ 931-1694/1117) -32.5 W32-W43 N08-N16
      7 1462  941 (N03W27 L=150)   452    166        1 (1414/ 894-1521/1017) -26.6 W22-W31 N01-N05
      8 1873  952 (N16W60 L=183)   546     90        2 (1850/ 921-1894/ 985) -63.4 W57-W63 N14-N18
      9  641  903 (S21E19 L=103)   459    219        0 ( 543/ 863- 715/ 958)  45.0 E25-E15 S26-S16
     10  891  817 (S20E02 L=121)   297     96        1 ( 828/ 793- 936/ 839)  38.7 E05-W00 S21-S17
     11  452  790 (S33E31 L= 92)   246     68        0 ( 429/ 765- 478/ 817)  59.0 E33-E29 S35-S30
     12 1262  611 (S21W26 L=149)   618    150        0 (1190/ 571-1345/ 656)  53.1 W23-W29 S25-S16
     13 1014  618 (S28W11 L=134)   684    137        1 ( 936/ 598-1083/ 655)  32.0 W07-W15 S31-S24
     14 1592  527 (S15W52 L=175)  1136    264        3 (1478/ 469-1679/ 595)  50.2 W45-W59 S21-S09
     15 1606  452 (S19W58 L=181)   990    276        1 (1519/ 429-1693/ 477)  21.8 W50-W66 S22-S14
     16 1483  334 (S29W56 L=179)   665     99        3 (1438/ 310-1541/ 359)  39.8 W53-W60 S31-S26
     17 1245  240 (S44W48 L=171)   996    166        0 (1196/ 213-1306/ 275)  39.8 W43-W51 S46-S39
  

    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