Color  xc   yc  (LatLon LCarr)  Area Length Branches (  x0/y0     x1/y1  )  Tilt Location
      1 1397 1516 (N25W39 L=271)  1057    268        0 (1283/1490-1515/1553) -12.5 W30-W48 N21-N27
      2  816 1490 (N40E01 L=231)   664    175        1 ( 741/1463- 896/1516) -15.3 E07-W04 N36-N42
      3  605 1479 (N45E18 L=214)   245     50        2 ( 579/1462- 633/1501) -38.7 E20-E15 N43-N47
      4 1755 1470 (N11W67 L=299)   554     80        1 (1715/1453-1791/1487)  -9.5 W60-W72 N09-N11
      5 1497 1422 (N17W41 L=273)   304     55        1 (1474/1399-1522/1451) -90.0 W39-W42 N14-N18
      6  918 1236 (N22E01 L=231)   802    314        0 ( 777/1171-1031/1285) -40.2 E09-W04 N15-N26
      7  410 1103 (N25E39 L=193)  1760    397        4 ( 251/1014- 546/1228) -34.7 E55-E29 N17-N35
      8  387  867 (N12E45 L=188)   253    113        0 ( 354/ 838- 442/ 883)   7.1 E47-E39 N10-N12
      9 1543  725 (S24W25 L=257)   798    243        1 (1462/ 659-1637/ 821)  14.9 W18-W33 S26-S20
     10 1419  704 (S22W16 L=248)   805    204        1 (1378/ 611-1442/ 793)  59.0 W11-W18 S26-S16
     11  737  742 (S02E23 L=209)   339     79        1 ( 705/ 712- 758/ 782)  31.0 E26-E21 S04-N00
     12 1755  541 (S47W49 L=281)   855     94        1 (1714/ 521-1799/ 564) -19.4 W41-W58 S49-S42
     13 1589  416 (S50W28 L=260)   707     89        2 (1561/ 381-1616/ 452)  59.0 W24-W31 S52-S46
     14 1507  378 (S50W18 L=251)   451     40        1 (1484/ 352-1525/ 405)  78.7 W16-W20 S52-S46
     15  459  372 (S21E62 L=170)   542     57        1 ( 435/ 347- 490/ 393)  15.3 E67-E56 S22-S19
     16 1237  339 (S42E06 L=226)   557     55        1 (1208/ 311-1273/ 361) -21.8 E08-E02 S44-S39
     17 1336  269 (S53E01 L=231)   410     75        1 (1303/ 250-1374/ 290)  -9.5 E03-W03 S55-S51
  

    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