Color  xc   yc  (LatLon LCarr)  Area Length Branches (  x0/y0     x1/y1  )  Tilt Location
      1 1041 1513 (N32E05 L=287)   362     75        1 (1021/1483-1064/1547) -76.0 E06-E02 N30-N34
      2 1287 1206 (N15W15 L=307)   449    142        1 (1242/1155-1316/1277) -50.2 W11-W16 N11-N18
      3 1641 1126 (N13W41 L=333)  1134    270        3 (1578/1005-1699/1235) -55.3 W35-W45 N06-N19
      4 1854 1054 (N11W63 L=354)  2254    306        3 (1816/ 904-1877/1192)  77.5 W58-W66 N01-N19
      5  143 1008 (S10E71 L=221)   956     94        0 ( 100/ 985- 187/1042)  -3.0 E82-E63 S10-S08
      6  789  890 (S09E13 L=279)   762    185        1 ( 749/ 810- 836/ 960)  63.4 E15-E10 S14-S04
      7  314  862 (S17E49 L=243)  1270    186        3 ( 247/ 780- 357/ 927) -66.0 E56-E44 S21-S12
      8  196  804 (S22E65 L=227)  1967    149        3 ( 163/ 736- 242/ 883) -84.3 E70-E58 S26-S16
      9  283  746 (S25E54 L=238)  1100    146        3 ( 260/ 676- 307/ 815) -83.7 E57-E51 S29-S20
     10 1530  582 (S21W41 L=333)   566    170        1 (1461/ 550-1574/ 628)  26.6 W35-W43 S23-S17
     11 1244  426 (S35W24 L=316)   682    172        0 (1231/ 356-1261/ 525) -69.0 W22-W28 S40-S27
  

    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