Color  xc   yc  (LatLon LCarr)  Area Length Branches (  x0/y0     x1/y1  )  Tilt Location
      1 1221 1528 (N27W27 L=313)   562    104        2 (1181/1497-1280/1548) -15.9 W24-W31 N25-N29
      2 1606 1529 (N15W53 L=339)   933    123        1 (1540/1510-1676/1548) -17.1 W47-W60 N13-N18
      3  249 1335 (N43E54 L=232)   649     96        1 ( 207/1303- 298/1360)   0.0 E60-E46 N41-N44
      4 1286 1281 (N12W22 L=307)  1236    434        0 (1126/1179-1454/1339) -41.2 W12-W28 N01-N18
      5 1748 1278 (S02W52 L=338)   778    132        0 (1686/1241-1792/1303) -18.4 W46-W56 S04-N00
      6 1851 1178 (S10W59 L=345)  1396    217        3 (1824/1082-1902/1254)  41.6 W54-W68 S15-S06
      7  799 1111 (N16E10 L=275)   499    170        1 ( 771/1037- 846/1184) -90.0 E11-E09 N10-N20
      8  313 1003 (N21E47 L=239)   332     85        0 ( 288/ 976- 348/1042) -51.3 E48-E43 N18-N23
      9 1702  779 (S28W37 L=322)  2355    517        0 (1565/ 555-1826/ 946)  30.7 W22-W50 S39-S21
     10  698  892 (N06E22 L=264)   626    118        3 ( 642/ 865- 767/ 928) -29.7 E24-E17 N03-N08
     11  271  886 (N15E53 L=233)   310     50        0 ( 247/ 870- 300/ 900) -18.4 E55-E49 N14-N16
     12 1302  642 (S24W06 L=291)  1041    332        3 (1206/ 532-1406/ 758)  25.2 E03-W14 S27-S19
     13  704  583 (S11E31 L=255)  1504    371        2 ( 634/ 515- 779/ 687) -90.0 E33-E25 S16-S04
     14  389  446 (S13E61 L=224)  5050    450       16 ( 259/ 309- 522/ 560) -76.0 E73-E49 S23-S03
     15  543  509 (S12E44 L=241)   369     76        0 ( 512/ 470- 572/ 536) -90.0 E46-E43 S14-S09
     16 1429  396 (S45W08 L=294)   319     83        0 (1403/ 359-1450/ 439)  45.0 W05-W10 S47-S42
     17 1005  315 (S36E25 L=261)  1019    131        2 ( 952/ 274-1065/ 359)  10.3 E29-E18 S39-S33
     18  888  151 (S46E49 L=237)   559     50        0 ( 863/ 123- 909/ 181)  45.0 E54-E44 S49-S42
  

    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