Color  xc   yc  (LatLon LCarr)  Area Length Branches (  x0/y0     x1/y1  )  Tilt Location
      1 1187 1867 (N59E25 L=193)   642     69        1 (1158/1848-1224/1886)  63.4 E29-E22 N55-N62
      2 1321 1787 (N56E06 L=213)   434     79        0 (1293/1762-1349/1806) -14.0 E09-E01 N53-N58
      3 1030 1568 (N28E16 L=203)  1096    305        0 ( 904/1494-1153/1640)  -5.4 E26-E05 N26-N29
      4  386 1594 (N14E64 L=155)  1060    125        4 ( 344/1536- 438/1624)  35.5 E68-E59 N09-N16
      5 1490 1416 (N32W19 L=238)  1345    216        3 (1428/1331-1547/1510) -33.7 W12-W24 N27-N36
      6  970 1419 (N17E14 L=205)   208     92        1 ( 925/1393-1007/1441)  21.8 E16-E11 N15-N18
      7  359 1345 (S02E52 L=167)   333     32        0 ( 344/1321- 375/1373) -36.9 E54-E50 S03-N00
      8  348 1218 (S09E49 L=170)   892    139        1 ( 304/1158- 392/1262)  54.5 E51-E45 S13-S06
      9 1460 1113 (N12W23 L=242)   650    158        0 (1441/1041-1475/1201) -71.6 W21-W24 N08-N17
     10  222 1068 (S21E58 L=160)   442     62        0 ( 184/1054- 255/1083)  12.5 E63-E54 S22-S19
     11 1659  966 (N10W41 L=259)   288     95        0 (1644/ 919-1687/1005)  90.0 W39-W42 N07-N12
     12  269  812 (S36E52 L=167)   763    174        0 ( 203/ 758- 330/ 870)  29.1 E62-E44 S40-S30
     13 1168  705 (S18W17 L=236)   374     75        1 (1139/ 679-1211/ 737) -11.3 W15-W20 S19-S16
     14  324  645 (S46E47 L=172)   758    181        1 ( 256/ 613- 395/ 683)  15.3 E59-E37 S49-S42
     15 1402  510 (S22W40 L=259)   356    100        1 (1371/ 471-1440/ 555) -80.5 W38-W41 S24-S18
     16 1343  503 (S25W37 L=255)   294    119        0 (1299/ 465-1383/ 540)  71.6 W35-W37 S27-S21
  

    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