Color  xc   yc  (LatLon LCarr)  Area Length Branches (  x0/y0     x1/y1  )  Tilt Location
      1  584 1483 (N08E39 L=186)   405     86        1 ( 546/1459- 629/1509)  45.0 E41-E37 N06-N10
      2  455 1465 (N04E48 L=177)   500     91        1 ( 427/1445- 497/1490)   0.0 E50-E44 N02-N06
      3 1170 1419 (N19E02 L=223)   317     68        1 (1141/1393-1202/1442)  68.2 E03-E01 N16-N21
      4  270 1345 (S06E60 L=165)   494    116        2 ( 245/1313- 287/1386) -63.4 E63-E57 S07-S03
      5  213 1253 (S12E62 L=163)  1519    321        4 ( 183/1189- 249/1315) -49.4 E68-E58 S15-S08
      6  268 1177 (S16E54 L=171)   454     87        0 ( 232/1152- 292/1209)  53.1 E57-E51 S17-S13
      7 1569 1150 (N15W28 L=254)   225     56        0 (1550/1117-1587/1180) -26.6 W26-W30 N13-N16
      8  894 1101 (S06E09 L=216)   402    161        1 ( 812/1064- 962/1118)   5.7 E14-E04 S08-S04
      9  621  875 (S26E20 L=205)   219     96        0 ( 578/ 858- 668/ 889)  38.7 E23-E17 S27-S23
     10  408  833 (S33E36 L=189)   521    147        1 ( 346/ 816- 472/ 852)  10.3 E41-E30 S34-S31
     11  964  683 (S28W06 L=232)   318    122        1 ( 907/ 665-1025/ 710)   7.1 W02-W10 S29-S26
     12 1382  648 (S17W32 L=257)  1455    503        0 (1160/ 612-1593/ 685)  20.8 W18-W47 S22-S11
     13 1688  543 (S12W60 L=285)   378     91        1 (1654/ 520-1720/ 561)  26.6 W55-W62 S13-S10
     14 1586  476 (S19W55 L=280)   872    293        0 (1530/ 408-1626/ 543)  -7.6 W46-W62 S21-S16
     15 1422  490 (S24W42 L=267)   747    188        1 (1329/ 469-1499/ 516)  28.6 W35-W46 S26-S20
     16 1292  350 (S36W42 L=268)   516    125        1 (1230/ 335-1349/ 370)  18.4 W37-W46 S38-S34
  

    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