Color  xc   yc  (LatLon LCarr)  Area Length Branches (  x0/y0     x1/y1  )  Tilt Location
      1 1369 1495 (N28W18 L=318)   256     75        0 (1340/1455-1389/1518) -63.4 W15-W19 N25-N29
      2  629 1382 (N12E29 L=272)  1155    306        2 ( 552/1295- 709/1508) -42.7 E36-E22 N07-N20
      3  399 1380 (N11E46 L=255)   280     55        2 ( 380/1348- 419/1411) -76.0 E46-E44 N08-N12
      4 1798 1155 (N12W51 L=352)  1983    406        6 (1707/1048-1904/1221) -10.8 W43-W64 N07-N16
      5 1395 1111 (N04W21 L=322)   323     85        1 (1362/1081-1432/1138)  36.9 W19-W23 N01-N05
      6 1337 1021 (S02W19 L=320)   267    109        1 (1299/ 991-1389/1068) -15.9 W15-W22 S03-N00
      7  163  928 (S17E63 L=238)  1445    236        4 ( 114/ 843- 211/1005) -48.0 E69-E57 S22-S12
      8  658  895 (S17E22 L=279)   208     65        0 ( 622/ 883- 694/ 905)  14.0 E23-E19 S17-S15
      9 1711  780 (S11W49 L=350)   448     91        1 (1681/ 746-1744/ 810) -26.6 W46-W53 S12-S09
     10  205  761 (S27E59 L=242)   419    138        1 ( 177/ 741- 230/ 789)  23.2 E62-E55 S28-S25
     11 1342  616 (S26W26 L=327)  1622    397        1 (1178/ 539-1535/ 748)  38.7 W15-W36 S32-S15
     12 1041  514 (S37W08 L=309)   396     79        1 (1000/ 500-1089/ 539)   8.1 W04-W11 S38-S35
  

    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