Color  xc   yc  (LatLon LCarr)  Area Length Branches (  x0/y0     x1/y1  )  Tilt Location
      1 1613 1668 (N51W35 L=195)   886    165        1 (1567/1639-1643/1700)  26.6 W29-W40 N47-N55
      2  877 1402 (N12E16 L=144)   763    247        0 ( 858/1293- 903/1520) -61.7 E19-E12 N06-N19
      3  235 1362 (N01E63 L= 97)   635    131        2 ( 208/1333- 263/1397) -18.4 E68-E59 N00-N03
      4 1094 1340 (N13E02 L=158)   208     83        0 (1054/1300-1132/1362)  59.0 E03-W00 N09-N14
      5 1145 1246 (N08W03 L=163)   516    162        1 (1107/1177-1175/1335) -84.3 W01-W04 N03-N13
      6 1563 1297 (N20W28 L=188)   277     75        0 (1550/1259-1576/1335)  80.5 W27-W29 N17-N23
      7  497 1309 (N01E38 L=122)   236     32        0 ( 472/1291- 515/1325) -18.4 E39-E36 N00-N02
      8  351 1196 (S08E47 L=113)   330     26        1 ( 326/1178- 372/1221) -14.0 E48-E44 S08-S05
      9 1544  698 (S14W40 L=199)  1105    355        2 (1418/ 594-1651/ 825)  66.0 W34-W43 S22-S04
     10  970  721 (S25W03 L=163)   912    335        1 ( 776/ 666-1101/ 776)  37.4 E08-W09 S32-S19
     11 1244  603 (S27W23 L=183)   281     61        0 (1212/ 588-1278/ 621)  26.6 W21-W25 S28-S25
     12 1465  329 (S37W52 L=212)  2401    465        2 (1281/ 302-1634/ 367)  25.6 W39-W64 S43-S30
  

    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