Color  xc   yc  (LatLon LCarr)  Area Length Branches (  x0/y0     x1/y1  )  Tilt Location
      1  873 1805 (N49E36 L=255)   512     81        0 ( 844/1777- 914/1823)  59.0 E38-E32 N46-N51
      2  654 1739 (N39E49 L=241)   646    103        0 ( 634/1698- 677/1790) -24.8 E56-E43 N35-N42
      3 1570 1556 (N44W35 L=325)   565     90        2 (1525/1537-1615/1577)   0.0 W30-W40 N42-N45
      4 1411 1439 (N32W19 L=309)   273     60        0 (1386/1416-1441/1460) -26.6 W17-W21 N31-N33
      5 1103 1355 (N21E01 L=289)   360    126        0 (1061/1305-1155/1392) -29.7 E04-W03 N18-N22
      6  296 1205 (S02E54 L=237)   858    129        2 ( 245/1162- 342/1255) -16.7 E60-E48 S03-N00
      7 1788  644 (S10W66 L=356)   571     74        0 (1747/ 621-1827/ 664)  26.6 W61-W71 S12-S08
      8  957  560 (S30W04 L=294)   333     75        0 ( 915/ 546- 997/ 575)  21.8 W01-W06 S31-S28
      9 1011  380 (S42W13 L=304)   310     53        0 ( 983/ 363-1032/ 402)  45.0 W11-W15 S43-S40
     10 1527  332 (S35W61 L=351)   575     62        2 (1505/ 310-1546/ 364) -76.0 W58-W66 S36-S32
     11 1576  336 (S34W67 L=358)   603     42        1 (1555/ 318-1599/ 364) -15.3 W62-W73 S35-S31
  

    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