Color  xc   yc  (LatLon LCarr)  Area Length Branches (  x0/y0     x1/y1  )  Tilt Location
      1  650 1144 (N20E22 L=197)   575    197        1 ( 569/1111- 737/1174) -14.0 E27-E15 N17-N22
      2 1567  932 (S09W32 L=251)   660    124        3 (1523/ 882-1604/ 989) -76.0 W29-W34 S12-S04
      3  762  693 (S09E22 L=197)   481    106        2 ( 725/ 657- 789/ 737)  68.2 E24-E19 S11-S06
      4  352  694 (S04E52 L=166)   392     64        1 ( 323/ 664- 371/ 724) -76.0 E54-E50 S06-S02
      5  296  603 (S10E62 L=156)   934    134        1 ( 265/ 553- 335/ 667) -82.9 E64-E59 S13-S05
      6  502  541 (S16E45 L=174)  1254    329        4 ( 439/ 476- 555/ 652)  61.4 E50-E39 S19-S08
      7  339  627 (S08E56 L=162)   510     65        1 ( 315/ 597- 369/ 649)  26.6 E59-E52 S10-S06
      8 1659  560 (S36W40 L=258)  1389    180        2 (1569/ 524-1729/ 607)   0.0 W30-W47 S38-S34
      9 1390  512 (S33W16 L=234)   689    152        0 (1327/ 483-1446/ 557)  11.3 W10-W20 S34-S30
     10 1021  454 (S29E11 L=207)  1278    319        1 ( 915/ 363-1126/ 553)  27.8 E21-E02 S34-S24
     11  684  479 (S22E33 L=185)   407     74        2 ( 656/ 447- 723/ 509)   9.5 E35-E29 S24-S19
     12  549  333 (S31E54 L=164)  2091    253        2 ( 517/ 239- 601/ 428)  26.6 E75-E43 S40-S24
  

    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