Color  xc   yc  (LatLon LCarr)  Area Length Branches (  x0/y0     x1/y1  )  Tilt Location
      1 1167 1647 (N45W08 L=263)   560     89        4 (1130/1609-1211/1673)  15.9 W04-W12 N41-N46
      2 1441 1606 (N43W32 L=287)  1640    281        2 (1301/1582-1555/1632)  -2.4 W19-W44 N40-N45
      3  405 1422 (N23E49 L=206)   461     77        1 ( 370/1400- 437/1453)  31.0 E51-E46 N21-N25
      4 1384 1134 (N11W22 L=277)  1954    406        5 (1321/ 945-1426/1320) -77.2 W18-W25 N00-N22
      5  650 1076 (N04E24 L=231)  1536    444        2 ( 600/ 897- 711/1264) -79.2 E28-E19 S07-N14
      6  830 1153 (N09E13 L=242)   374    120        2 ( 809/1094- 859/1203) -74.1 E14-E10 N05-N12
      7  165 1124 (N03E67 L=187)   631    108        0 ( 145/1077- 185/1186)  63.4 E69-E65 N00-N06
      8  288  904 (S09E51 L=203)   815    160        0 ( 266/ 832- 317/ 980) -77.5 E53-E48 S13-S04
      9  943  877 (S07E04 L=250)   657    211        0 ( 915/ 780- 984/ 965)  79.7 E05-E02 S12-S01
     10 1625  723 (S14W43 L=297)   328    140        0 (1585/ 695-1680/ 751)  -7.1 W39-W47 S15-S12
     11  872  656 (S21E08 L=246)   580    136        0 ( 825/ 613- 936/ 694)  32.0 E11-E03 S24-S18
     12 1536  633 (S20W38 L=292)   275     74        0 (1508/ 619-1578/ 654)  -8.1 W34-W41 S20-S18
     13  846  532 (S30E10 L=245)   411     60        0 ( 829/ 497- 869/ 577) -80.5 E11-E08 S32-S26
     14 1393  300 (S46W40 L=295)   457     63        1 (1359/ 285-1427/ 319)   8.1 W36-W43 S47-S43
     15  662  259 (S55E35 L=220)  1099     93        2 ( 620/ 228- 700/ 292) -20.6 E40-E28 S59-S51
     16  737  264 (S54E24 L=230)   517     33        1 ( 714/ 242- 760/ 287)  26.6 E28-E20 S56-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